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Optical nanofibers provide a rich platform for exploring atomic and optical phenomena even 
when they support only a single spatial mode. Nanofibers supporting higher-order modes provide 
additional degrees of freedom to enable complex evanescent field profiles for interaction with the 
surrounding medium, but local control of these profiles requires nondestructive evaluation of the 
propagating fields. Here, we use Rayleigh scattering for rapid measurement of the propagation 
of light in few-mode optical nanofibers. Imaging the Rayleigh scattered light provides direct vi¬ 
sualization of the spatial evolution of propagating fields throughout the entire fiber, including the 
transition from core-cladding guidance to cladding-air guidance. We resolve the interference between 
higher-order modes to determine local beat lengths and modal content along the fiber, and show 
that the modal superposition in the waist can be systematically controlled by adjusting the input 
superposition. With this diagnostic we can measure variations in the radius of the fiber waist to 
below 3 nm in situ using purely optical means. This nondestructive technique also provides useful 
insight into light propagation in optical nanofibers. 


I. INTRODUCTION 

Recently, optical nanofibers (ONE) have received sig¬ 
nificant attention by realizing efficient interactions with 
atoms [MU and chemical analytes [l^ - [T3| . Knowledge 
of light propagation in ONFs is particularly important 
for atomic systems, for which the strong evanescent field 
enables both tightly trapped atoms near the nanofiber 
and efficient atom-light coupling for interrogation. For 
a single-mode ONF that supports only the fundamental 
optical mode, the evanescent field distribution is well- 
defined and controlled; these ONFs have produced im¬ 
portant results in atom-light interactions [lH6|, but have 
been limited to only a few optical potential configura¬ 
tions. 

To increase control over the evanescent field of the 
ONF and the complexity of the potential landscape, re¬ 
cent studies have explored the use of higher-order vec¬ 
tor modes [l5l - ll7ll. which have spatially-varying polar¬ 
ization profiles [l8l Although the ability to cre¬ 

ate, control, and detect complex optical mode distri¬ 
butions [H [H [iii3 at the input and output of a 
fiber is well-established, the challenge remains to ana¬ 
lyze the propagation of light along the ONF taper and to 
measure the local optical intensities on the ONF waist. 
In a single-mode ONF, the polarization dependence of 
Rayleigh scattering can discriminate between the two 
degenerate, quasi-linearly-polarized, fundamental modes 
along the fiber [lOl . l23l - l^ , but such discrimination is not 
possible for higher-order modes (HOMs) without high 
transverse spatial resolution Additionally, because 
ONFs are typically drawn to radii below the resolution 


limits of optical imaging, the ONF geometry is difficult to 
quantify in situ. Propagation dynamics through the fiber 
can be inferred from analysis of the transmitted light dur¬ 
ing the drawing process [ilH , but local visualization 
of the modes has not yet been achieved. 

In this paper, we image Rayleigh scattered (RS) light 
to analyze light propagation in ONFs supporting the first 
six optical modes. By combining high-spatial-resolution 
imaging with the polarization dependence of RS, we mea¬ 
sure and control the propagation along the entire ONF 
length. The intermodal beat lengths, which vary along 
these axially-tailored fibers, are quantified over the ONF 
and are sensitive to changes in the fiber radius to within 
3 nm. Finally, we show that by controlling the input 
modal superposition, we can systematically adjust the 
modal content on the ONF waist. 

This paper is organized as follows. Section HIl describes 
the dispersion relations in ONFs and properties of the 
modes of interest. Section m details the experimen¬ 
tal setup used to fabricate and analyze propagation in 
the ONFs. In Sec. M we present the key experimen¬ 
tal results. We present RS images over an entire ONF, 
including the transition from core-cladding guidance to 
cladding-air guidance. We demonstrate modal identifica¬ 
tion and control of the mode distribution. We also dis¬ 
cuss mode identification by further polarization discrim¬ 
ination of the RS light into longitudinal and transverse 
components. Section 5 concludes the paper. 
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FIG. 1. Geometry of an optical nanofiber. The camera view¬ 
ing direction is along y, and light propagates from left to right 
along z. E is the scattered field. The tapered core is negligible 
in the fiber waist. 


II. THEORY 

A. Nanofiber geometry and modes 

An ONF consists of three sections: the unmodified 
fiber, the taper, and the fiber waist (See Fig. [T]), with 
a symmetric taper on the output. Initially we launch 
light into the core of the unmodified fiber on the left. 
The taper makes an angle O with the fiber axis and adia- 
batically links the unmodified fiber to the fiber waist [2^ . 
Near the waist, the taper angle is reduced exponentially 
to connect the taper to the waist. The waist has a sub¬ 
micrometer radius, and length, that is typically 
several millimeters long. 

The number of guided fiber modes depends on the ra¬ 
dius, a, the wavelength of input light, A, and the refrac¬ 
tive indices of the core and cladding, ni and n 2 , through 
the dimensionless lA-number: 


V = ak^^nl -nl, (1) 

where k = 27 r/A is the free-space wave vector. The al¬ 
lowed modes and their effective indices, Ueff are shown 
in Fig. H a) for an ONF with 1^ < 8 (a < 0.9 /im). In 
a tapered fiber, V changes with propagation distance. 
Previous works have examined taper geometries that al¬ 
low adiabat ic p ropagation to minimize interfamily mode¬ 
mixing [I^, [ 23 . In this work, the ONFs are drawn from 
Fibercore SM1500 optical fiber, which has V = 3.8 for 
our operating wavelength A = 795 nm, and initially sup¬ 
ports the LPqi, LPii, and LP 21 mode families. The LPu 
family is composed of the TEqi, TMqi, HE 21 and HE 21 
vector modes, and their transverse field distributions are 
shown in Fig. HKb). 

Each mode has an associated propagation constant 
/3 = Ueffk. Two propagating modes interfere with a 
beat length, Zf, defined by 


Zb{a) = 


2tt 


A 


/ 3 i(a) -/32(a) neff,i(a) - neff,2(a) 


( 2 ) 


As we will show, the evolution of Zb along the length of 
the fiber due to the tapered geometry is readily detected 
by imaging the RS light, and can be used to measure 
local fiber dimensions. 
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FIG. 2. The effective index of refraction for the lowest lying 
modes in a nanofiber versus V-number (bottom axis) and ra¬ 
dius (top axis), b) The vector profile inside a 450 nm radius 
fiber for the transverse components of the TEqi, TMqi and 
HE 21 modes respectively. The HE 21 mode is not shown, but 
is orthogonal at all points with the HE^i mode, c) HE 21 and 
TMoi modes interfering d) HE 21 and TEqi modes interfer¬ 
ing. These images are representative of the transverse fields 
on the ONF waist. The solid black circle marks the cladding 
air interface at a radius of 450 nm. The refractive indices are 
chosen to be 1.45 and 1 for the cladding and air respectively. 


The propagation constants in Fig. [2] have also been 
plotted as a function of a (top axis). From Eqn. [21 we 
can calculate the beat lengths for the modes of interest 
as a function of local fiber radius. For an ONF with « 
400 nm, Zb varies from a few microns for HOMs beating 
with the fundamental mode, to tens of microns for HOMs 
beating together. The TEqi and TMqi mode pair and 
the HE 21 and HE 21 mode pair do not interfere because 
their polarizations are orthogonal everywhere. 


B. Rayleigh scattering 

Standard optical fibers have birefringence due to 
strain, microbends, or deviations from cylindrical sym¬ 
metry that can convert a pure input mode to a superposi¬ 
tion of modes on the fiber waist. This necessitates a non¬ 
destructive, in situ measurement that locally probes the 
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intensity and polarization. Using Rayleigh scattering, we 
can extract information about the propagation and mode 
distribution along the ONF. The observed RS power is 
proportional sin^d, where 9 is the angle between the local 
field, E, and the camera viewing direction, y. When the 
mode-field-diameter is large enough, imaged RS shows 
the spatial structure of the propagating light and 
in Sec. IIVIIV Al this will be used to observe the transi¬ 
tion from core-cladding guidance to cladding-air guid¬ 
ance. On the waist of the ONF, which typically has 
ttu) < 400 nm, transverse mode structure cannot clearly 
be resolved. However, since Zb is much smaller than the 
length of the waist, the local superposition on the waist 
is determined straightforwardly through Eq. [5] 


III. EXPERIMENTAL SETUP 
A. Fiber pulling 

We use the flame brush technique to produce the 
ONFs [ISlIIl. For the full details of the pulling appara¬ 
tus see Refs, [l^ [H, [2^. It consists of two computer- 
controlled high-precision motors and an oxyhydrogen 
flame in a stoichiometric mixture. We employ an algo¬ 
rithm [28l - l3Tj | to produce ONFs with linear tapers. Fol¬ 
lowing the protocols outlined in Ref. [1^ we can produce 
low-loss ONFs when controllably launching the funda¬ 
mental mode [ 1 ^, [ 2 ^ or the next family of higher-order 
modes [isj ]. Fiber cleanliness prior to tapering is crucial 
to achieving high transmission [ 2 ^. 


B. Mode preparation and imaging setup 

Figure |3] shows the experimental setup. It has three 
components: the ONF fiber launch and vector beam gen¬ 
eration setup [H, the RS imaging setup, and the 
ONF output imaging setup. The HOMs coupled into 
the ONF are derived from a fiber-based mode converter, 
described in detail in Ref. M- Briefly, we begin with 
spatially-filtered, polarized light from a 795 nm diode 
laser. This beam is collimated and passed through a 
phase plate that imparts a tt phase shift onto half of the 
beam, producing a two-lobed profile at the entrance of a 
vector beam generator. 

The vector beam generator uses fiber that supports 
only the HEn, Ti?oi, TMgi, and 77 £'21 modes at 795 
nm (Corning HI1060). A polarization controller (PC) 
mixes these modes through strain-induced birefringence 
and can produce either pure modes or any superposition 
of them at the output |^. Additional control over the 
mode superposition is done with bulk polarization optics 
(e.g. A/2 waveplate in Fig.|31) at the input of the ONF. A 
CCD at the output of the ONF is used to aid alignment 
into the ONF. The overall fiber length is 20-30 cm, and 



FIG. 3. The experimental setup to measure the Rayleigh 
scattering from an ONF. 

it is kept straight and tensioned slightly to reduce drifts. 
We do not observe fluctuations over periods of several 
hours. 

After we fabricate an ONF we use the motor stages 
(two high precision XML 210 motors see Fig. [3]) to trans¬ 
late the fiber relative to the RS imaging system to im¬ 
age the entire ONF. RS light emitted orthogonal to the 
propagation direction is relayed to an EMCCD camera 
(Andor Luca) through a Mitutoyo M Plan APO NIR 
infinity-corrected objective (lOx or 50x). The imaging 
relay contains a 795 nm bandpass filter and a calcite 
walk-off polarizer near the CCD image plane to sepa¬ 
rate the longitudinal- and transverse-polarized RS light. 
When using the calcite polarizer, the spatial resolution 
of the imaging system is limited to 2-3 p,m. Because the 
beat length between the LPn and LPqi families is also 
2-3 p,m, this interference is not observed on the ONE 
waist. 


IV. RESULTS 

A. Direct imaging of optical propagation 

In a typical fiber, light is confined to the fiber core 
by the lower-index cladding, and the air surrounding 
the cladding has a negligible effect on the propagation. 
On the ONF waist, however, the air itself becomes the 
cladding and the original fiber core is negligible. The 
transition between these two regimes has previously only 
been indirectly studied by analyzing the transmitted 
power [ 3 , • In this section, we describe the di¬ 

rect imaging of the propagation evolution throughout the 
ONF. 

Using the experimental setup described in Sec. Ell we 
launch a HOM into an ONF. To observe the propagation 
over the entire 10-cm fiber length, we use the lOx objec¬ 
tive, which has a field-of-view of Ril mm. Figure HJa) 
shows a RS image produced after stitching together the 
series of images covering an ONF with = 0.75 mrad and 
Ou) = 1.5 /im. In this figure, light is propagating from left 
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FIG. 4. High-resolution RS images of an ONF. Light is prop¬ 
agating from left to right, (a) Propagation through the entire 
ONF, beginning and ending in unmodihed hber. = 0.75 
mrad, = 1.5pm. b) Magnihcation of the input taper 
showing the transition from core to cladding near a radius 
of 13.5 pm. c) Simulation of the propagation region shown 
in (b). Colors in the images have been adjusted to highlight 
structure in the taper. 


to right. The light that is initially in the core is ejected 
into the cladding when the fiber cladding has a diameter 
of 27 pm. This is consistent with the value obtained in¬ 
directly in Ref. [l^. On the output side, the field couples 
back into the core at the same radius, but the RS power 
is reduced due to light that escaped from the ONF on 
the input taper. This light is lost because the ONF waist 
only supports the LPqi and LPn mode families, whereas 
the unmodified fiber initially supported even higher order 
mode families at the fiber launch. Additionally, higher- 
order mode families could be excited due to nonadiabatic 
propagation, as we have previously found that adiabatic- 
ity for the LPn family requires O < 0.4 mrad M- The 
color scaling in the image has been adjusted to highlight 
details in the taper and does not indicate saturation. 

Figure. IHb) is a magnification of Fig. IH^a) over 4 
mm of propagation distance in the critical region near 
a = 13.5 pm, where the ROMs transition from the core 
to the cladding. For these images we used a 50x objec¬ 
tive (NA = 0.65) to collect the transverse polarization 
component of the RS light. The core mode shows a null 


in the center of the scattering, which is an expected con¬ 
sequence of HOM Rayleigh scattering as superpo¬ 
sitions of ROMs can have a two-lobed intensity pattern. 
It is also possible that the field is a pure HE 21 or TMqi 
mode, which would also appear with a null because of the 
polarization filtering of RS . After the core-cladding 
transition, the oscillations observed in Figs. HJa) and (b) 
increase in spatial frequency as the fiber tapers down to 
the ONF waist. As will be discussed in Sec. IIVIIVDI 
the total RS power increases toward the waist because of 
surface scattering effects. 

We can analyze the length scales associated with the 
transition from core-cladding guidance to cladding-air 
guidance using the 1/e decay of the scattered power in 
the core. This decay occurs over a propagation distance 
of ~700 pm. For a fiber with R = 0.75 mrad this cor¬ 
responds to a cladding radius change of 500 nm. The 
specifications of SM1500 hber indicate that the hber has 
an initial core radius of 1.8 pm and cladding radius of 
25 pm. Tapering does not affect the core-to-cladding ra¬ 
tio, so this transition occurs over a change in core radius 
of ~40 nm. The oscillation pattern in the transition re¬ 
gion is consistent with the length scales associated with 
the LPii family beating with the next symmetric family 
LP 21 . In this region, the relevant F-number governing 
the propagation increases rapidly from V 1.9 before 
the transition to F r; 115 after the transition. Recause 
the number of allowed modes scales with F^, this is the 
region where intermodal transitions are most likely to oc¬ 
cur and the adiabaticity criteria are most stringent (^ . 

Figure. [D)c) is a simulation using commerical waveg¬ 
uide propagation software (FIMMPROP/FIMMWAVE 
by Photon Design, Ltd.) of the local intensity, integrated 
along the viewing direction, over the region in Fig. H^b). 
This calculated image is for a TMqi mode, and does not 
include the polarization effects of RS, but the agreement 
between experiment and simulation is clear, including the 
location at which the light leaves the core and the chirped 
oscillation frequencies observed in the cladding. In this 
region, the number of modes is large enough that it is pos¬ 
sible to think in a pure ray optics picture, for insight, to 
illustrate the spatial oscillations of light and their change 
of frequency. 


B. Mode identification 

We can spatially resolve transverse modal information 
from RS when the light is confined to the core and when 
the ONF is larger than 3 pm, but direct imaging on the 
ONF waist is difficult, where Qw ~ 400 nm. Rowever, the 
oscillations observed in Fig. S] are due to superpositions 
of modes that can also be identified by spatial frequency 
analysis of RS. Such analysis is used in standard, single¬ 
mode birefringent hbers [^. In the ONFs, however, sev¬ 
eral modes interfere with beat lengths that evolve along 
the fiber as the ONF radius, a, tapers. We determine 
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these Zb through Fourier transform of the RS power as 
a function of z. We use spectrogram analysis to 

obtain the inverse beat frequencies, l/zb- The spectro¬ 
gram shows all the spatial frequencies that evolve as a 
function of propagation distance (and therefore radius) 
in the fiber. 




FIG. 5. Calculated inverse beat frequency for = 360 nm 
with 17 = 1 mrad and = 5 mm as a function of a) propa¬ 
gation distance; and b) ONF radius. 

Figure [SKa) shows the calculated inverse beat lengths 
as a function of position for a fiber with a profile a^, = 360 
nm, 17 = 1 mrad, and = 5 mm using FIMMPROP. 
This is the shape we expect from the spectrogram for the 
HE 21 and TEqi (red) beating together and the HE 21 and 
TMqi (blue) interfering. As shown below, the large quan¬ 
titative and qualitative differences between these two 
curves enables us to identify the modes observed in the 
experiment. The spatial frequencies for both superposi¬ 
tions have also been plotted as functions of a in Fig.jSKb). 

Fig. [HKa) shows the power scattered as a function of z 
near the ONF waist for a fiber with 17 = 1 mrad and a 
desired = 390 nm. The spatial variation in oscillation 
frequency of the RS power is depicted in the spectrogram 
in Fig.lHKb), in which FFTs of the spatial domain data are 
calculated using 0.5 mm windows. The qualitative shape 
in Fig.lHKb) matches the shape of the HE 21 : TMqi curve 
from the calculation in Fig. [S^. The spatial frequency is 
zero at the crossing of the UeB curves for the HE 21 and 
TMqi modes near 430 nm, shown in Fig. [2J 

The steep dependences of 1/zb on a are shown in 
Fig. [5Kb) and allow a high-precision, local measure of 
variations in a. We extend the FFT window over the en¬ 
tire waist to increase the spatial frequency resolution and 
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FIG. 6. a) Detected (blue) and smoothed (red) transverse- 
polarized RS power along the ONF waist, b) A spectrogram 
of concatenated images displaying the HE 21 mode beating 
with the TMqi mode in arbitrary false color for the intensity. 
The dashed black line corresponds to values expected from a 
fiber profile calculated by the pulling software using a™ = 360 
nm. 


obtain a beat frequency 1/zf, = 24.50 ± 0.06 mm“^ with 
a full-width at half-maximum (FWHM) of 0.24 mm“^. 
Considering the ONF as a single dielectric core sur¬ 
rounded by an infinite cladding of air, we estimate from 
this beat frequency that = 360 ± 10 nm. With only 
the HE 21 '■ TMqi measurement, we cannot constrain 
with better than 10 nm uncertainty as it depends slightly 
on the exact index of the ONF. However, the measured 
width of 1/zb = 0.24 mm“^ limits the variation in 
to 0.7 nm over the entire 5 mm waist [see Fig. EKb)]. 
Detailed calculations show that a similar high-resolution 
measurement of the beat frequency between HE 21 and 
TEqi modes on the same fiber should constrain a^), and 
this will be considered elsewhere. The estimated value of 
Ou) = 360 nm is ~30 nm smaller than the design radius 
of 390 nm. Superimposed on the data in Fig. |5|is a beat 
frequency curve using a profile from the design software 
with a revised = 360 nm. 

Fig. 13(a) shows the RS profile for a different ONF sam¬ 
ple than the data displayed in Fig. [51 In this case, how¬ 
ever, the shape of the spectrogram curve indicates beat¬ 
ing between the HE 21 and TEqi modes. We follow a 
similar procedure for this ONF sample as the previous 
sample to extract a center beat frequency 1/Zb = 75.6 ± 
0.1 mm“^, with FWHM = 0.7 mm“^. Following a simi¬ 
lar procedure as above, we estimate = 370 ± 10 nm. 
At this radius, the beat frequency for the HE 21 ■ TEqi 
superposition has a weaker dependence on than the 
HE 21 : TMqi superposition and gives an estimated uni¬ 
formity of 3 nm over the 5 mm length of the waist. Beat 
frequency analysis provides a precise nondestructive ex¬ 
traction of the fiber profile without requiring high NA 






























6 


optics or irreversible coatings for imaging in a scanning 
electron microscope (SEM). 



FIG. 7. a) Detected (blue) and smoothed (red) transverse- 
polarized RS power along the ONF waist after summing over 
the transverse scattered power, b) A spectrogram of aligned 
concatenated images displaying the HE 21 mode beating with 
the TEoi mode in arbitrary false color for the intensity. The 
dashed black line corresponds to an overlapped exponential 
fiber profile. The overlap results in a hber radius of 370±10 
nm. Note this is for a different sample than the fiber displayed 
in Fig. [6l 


HWPs and polarizers, to launch known mode superposi¬ 
tions into the ONF and control the mode propagating on 
the ONF waist. 

Figure [2](b) displays the vector profile for the trans¬ 
verse components of a 450 nm radius fiber for the HOMs 
in the LPn family. In free space, conversion between 
any of the pure modes is readily accomplished with one 
or two HWPs. For example, by placing two HWPs with 
a relative angle of 45 degrees between their fast axes 
in front of a radially-polarized mode will produce an 
azimuthally-polarized mode at the fiber input. Likewise, 
the HE 21 polarization profile is formed when a single 
HWP is placed in the path of a TMqi profile. 

Similarly, when launching an equal superposition of the 
HE 21 and TMqi modes, which produces a two-lobed pro¬ 
file with polarization perpendicular to the null, we can 
completely convert to an equal superposition of E[E 2 i 
and TEqx using a HWP with an angle of 45 degrees 
with respect to vertical. Figs. EIc) and [2Kd) show the 
transverse profiles of the HE 21 : TMqi superposition 
and the HE 21 : TEqi superposition, respectively. Be¬ 
ginning with the HE 21 : TMqi distribution in Fig. He), 
it is straightforward to see that adding a HWP at angle 
a produces the HE 21 : TEqi superposition with relative 
power sin^2Q; and the HE 21 : TMqi superposition with 
relative power cos^2a. In a fiber, these superpositions 
may be modified due to strain or deviations from cylin¬ 
drical symmetry. 


The total scattered power increases by lOx in Fig. [5] 
and 20x in Fig.[7]as a decreases, reaching a maximum on 
the ONF waist. Because RS is directly proportional to 
intensity, the contribution from bulk RS should not in¬ 
crease. As a decreases, however, the field at the fiber sur¬ 
face increases, and it is likely that the observed increase 
in RS power comes from surface scattering. This scatter¬ 
ing appears as noise in the spectrograms and occurs at all 
spatial frequencies [see Fig.|6l[b) and Fig.[7l[b)]. The sur¬ 
face scattering also increases at radii where HOMs reach 
cutoff. This is apparent in Fig. where the signal 
drops sharply near z = 2 mm, at which point one of 
the LP 21 modes escapes the ONF. On the output side of 
the fiber, this RS signature is absent because this HOM 
was no longer present. Surface scattering will also be 
discussed in Sec. IIMIVDI 


C. Mode control 
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The ability to measure spatial frequencies and identify 
participating modes suggests that external control of the 
input mode distribution can be used to create a desired 
mode superposition on the waist. To demonstrate con¬ 
trol we launch and manipulate the manifold of HOMs in 
a systematic way. Following Ref. 0 we use the cylindri¬ 
cal vector beam generator and bulk optics, consisting of 


FIG. 8. Conversion from the HE 21 and TMqi modes inter¬ 
fering to HE 21 and TEqi. a) Each column corresponds to 
a FFT at a given HWP angle, b) The absolute value of the 
FFT along the rows in (a) for the HE 21 and TMqi beating to¬ 
gether (blue squares) and the HE 21 and TEqi modes beating 
together (red circles). Fits to the expected sin^2Q; are shown 
as solid lines. 
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Figure [8] shows the relative powers of these superposi¬ 
tions as a function of a. Beginning with the HE 21 : TMqi 
input superposition, we rotate a HWP in 2 degree incre¬ 
ments and collect the transverse RS over the ONF waist. 
Then we take the absolute value of the Fourier transform 
at each a and plot it as the column of an image, as shown 
in in Fig. |8l[a). Note that the longitudinal (z-polarized) 
RS cannot detect the HE^x '■ TEqx superposition, be¬ 
cause the TEqi mode is entirely transverse. 

The image shows two primary frequency bands, near 
1/zh = 75 mm”^ and l/z;, = 24 mm“^, close to the 
two values shown in Fig. [5] on the waist. The powers 
of these individual frequency bands are also plotted in 
Fig. [8jb), with the HE 21 : TMqi shown in blue and the 
HE21 ■ TEqi shown in red. As described above, the 
data fit well to a sin^2a dependence, indicated by the 
solid lines. At a = 0, most of the signal is due to the 
HE 21 : TMqi superposition. As a increases, the higher 
frequency band corresponding to the HE 21 : TMqi su¬ 
perposition increases and lower band corresponding to 
the HE 21 : TEqi superposition decreases. At a ~ 40 
degrees, the HE 21 : TEqi superposition is maximized on 
the ONF waist. Rotating the FIWP another 45 degrees 
to 85 degrees minimizes the HE2X : TEqi superposition. 
The two bands are out of phase with almost full con¬ 
version from one pair to the other. The small difference 
in amplitude of the two curves is due to the change in 
coupling into the ONF as the waveplate is rotated. 

This control has been done entirely with bulk optics 
and a straight ONF. For a fiber with bends, however, 
the transfer matrix between the input set of modes and 
the waist set of modes can have off-diagonal terms, re¬ 
quiring the ability to have arbitrary state generation at 
the input. For vector modes, bulk optics alone cannot 
generate all possible input superpositions. Polarization 
controllers relying on strain-induced birefringence, such 
as that used in the vector beam generator, can create any 
superposition of vector modes [20| . 


D. Surface scattering and mode cutoffs 

In Sec. IIMIV Cl we measured the modal superposition 
on the waist as a function of the input state by mea¬ 
suring spatial frequencies in the transverse-polarized RS 
data. However, in these tightly confining waveguides, the 
propagating field can also have a longitudinal component 
of similar strength to the transverse component. In this 
section, we analyze both the longitudinal and transverse 
components of the RS to enable more specific mode iden¬ 
tification. This discrimination is done by adding a calcite 
walkoff plate near the camera plane to enable simultane¬ 
ous measurement of both polarizations. We also discuss 
the effects of surface scattering, which provides a useful, 
indirect measure of the field intensity at the nanofiber 
surface. 


As can be seen in Fig. [2l the HE 21 modes reach cut¬ 
off at a larger radius than the other modes in the LPu 
family. By pulling a fiber to a radius just below the 
HE 21 cutoff, the only propagating modes are the TEqi 
and TMqi modes. For these ONFs and 795 nm light, 
the HE 21 cutoff is reached at = 325 nm. Fig.[S]shows 
the longitudinal, transverse, and total RS power from the 
ONF as a function of position along the fiber axis for O 
= 1 mrad, Ou, = 300 nm, and Luj= 1 cm. By pulling to 
this radius, a number of phenomena are observed. 

From 2 : 0 — 10 mm on the fiber axis, the total scat¬ 

tered power is constant. That this region, a taper section, 
has a steady scattering rate, is indicative of bulk scatter¬ 
ing. After this point, from z = 10 mm to z = 24 mm, 
however, there is a steady increase in the scattering rate 
as the fiber continues to taper, consistent with surface 
scattering effects as discussed earlier. The intensity of 
the mode at the fiber surface increases as the radius de¬ 
creases and surface scattering begins to dominate over 
bulk scattering. This increase in scattered power was 
also observed in Figs. O and[71 

The RS light abruptly drops at z ~ 25 mm, indicated 
by the dashed line at z = 24 mm in Fig. [HI This is the 
point that the HE 21 mode reaches cutoff at a = 325 nm, 
corresponding to a H ~ 2.8. The drop in RS power is 
most pronounced in the longitudinal polarization com¬ 
ponent in Fig. m a). In each panel of Fig. O the RS light 
is constant between z = 25 mm and z = 35 mm. This 
region, denoted Lu,, is the 10-mm-long ONF waist. The 
contribution of the HE 21 mode to the total scattered 
power is shown in Fig. He), indicated by the red dashed 
lines. These levels are taken at two symmetric places on 
opposite sides of the waist with equal radii. 

On the output side of the waist, the total RS power 
increases slightly. This is again a consequence of surface 
scattering. For this ONF, = 300 nm and is also very 
close to the cutoff for the TMqi and TEqi modes. At 
this radius, the remaining hollow modes have a signifi¬ 
cant evanescent field and low surface field. As the radius 
increases on the output, the evanescent field reduces and 
the surface held increases, indicated by the larger RS 
power near z = 37 mm. The surface held reaches a max¬ 
imum near a = 310 nm, beyond which the held becomes 
more strongly conhned inside the hber and the surface 
scattering reduces. The RS power continues to drop as 
more of the held is conhned to the bulk of the ONF. 

The polarization discrimination of Fig. combined 
with the cutoh of the HE 21 mode, helps identify the re¬ 
maining amount of the TMqi and TEqi modes. This 
is because the TEqi mode has no longitudinal held and 
it should not contribute to any of the power seen in 
Fig. Ha). Since the HE 21 modes are ejected from the 
waveguide, and scattering still exists in the transverse 
direction but not in the longitudinal direction, we know 
that the majority of the total RS power on the waist is 
due to the TEqi mode. We note that there is some resid¬ 
ual longitudinal component that could be caused by a 
small TMqi or HEu background that is not observable 
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FIG. 9. Launching a superposition of HOMs into a fiber with 
= 1 mrad, a™ = 300 nm, and Lw = 10 mm, marked with a 
thin, dashed black line. The three panels correspond to scat¬ 
tering collected with (a) longitudinal polarization, (b) trans¬ 
verse polarization, and (c) total scatter as a function of length 
respectively. The long, thick dashed black line denotes the 
HE 21 cutoff and from that point we observe only the propa¬ 
gation of the TEoi mode. The short dashed red lines show the 
power loss from the HE 21 mode ejecting from the ONF. The 
continuous black lines designate the ONF waist. The power 
scattered from 0 mm to 25 mm corresponds to the input taper 
and the scattered power from 35 mm to 50 mm corresponds 
to the output taper. For clarity, the detected RS signal (solid 
blue) is smoothed (solid red) in each panel. 


in the spectrogram. It could also be due to deviations 


from pure dipolar scattering [s^. Ref. notes that 
the longitudinal scattering was significantly larger than 
expected, even in unmodified fibers with negligible lon¬ 
gitudinal fields. These results indicate that pure TEqi 
modes can be identified by minimizing the longitudinally- 
polarized RS light. 


V. CONCLUSION 

We have shown that Rayleigh scattering is an excellent 
diagnostic of propagation phenomena in axially-tailored 
nanofibers. Rayleigh scattering enables direct observa¬ 
tion of the propagation throughout an optical nanofiber, 
allowing us to demonstrate control of the mode superpo¬ 
sitions on the nanofiber waist. By measuring the vari¬ 
ation of the beat lengths, we can determine the local 
fiber radius to within 3 nanometers using entirely op¬ 
tical means. The ability to measure and identify modes 
in high transmission, few-mode ONFs opens the possibil¬ 
ity of controlling HOM-based evanescent traps for atoms. 
This tool could also be extended to integrated waveguide 
approaches to atom traps . 
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